Strongly enhanced shot noise in chains of quantum dots 
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We study charge transport through a chain of quantum dots. The dots are fully coherent among 
each other and weakly coupled to metallic electrodes via the dots at the interface, thus modelling a 
molecular wire. If the non-local Coulomb interactions dominate over the inter-dot hopping we find 
strongly enhanced shot noise above the sequential tunneling threshold. The current is not enhanced 
in the region of enhanced noise, thus rendering the noise super-Poissonian. In contrast to earlier 
work this is achieved even in a fully symmetric system. The origin of this novel behavior lies in a 
competition of "slow" and "fast" transport channels that are formed due to the differing non-local 
wave functions and total spin of the states participating in transport. This strong enhancement may 
allow direct experimental detection of shot noise in a chain of lateral quantum dots. 

PACS numbers: 73.63.-b, 73.23.Hk, 72.70.+m 



Introduction- Shot noise in mesoscopic devices such 
as point contacts and quantum dots has been the sub- 
ject of intensive research within the last decade [jj. Re- 
cently the interest has been intensified both theoretic ally 
0H&SIS0I1I3 and experimentally [H [H E3 
when the importance of shot noise for quantum comput- 
ing, nanomechanical systems and molecular devices had 
been recognized. However, experiments on these systems 
are difficult to perform, since one needs to detect the 
shot noise over the background of l/f- noise caused by 
fluctuations in the physical environment and measure- 
ment equipment. So far even for lateral quantum dots in 
semiconductor heterostructures, with far reaching con- 
trol of system parameters and advanced low temperature 
setups, no direct shot noise measurements have been re- 
ported [13]. 

In an array of quantum dots with Coulomb interac- 
tions, both the current / and the shot noise S depend 
on details of the electronic spectrum of the many-body 
states of the coupled system, as well as the coupling 
strengths of these states to the electrodes 0, 0. Thus 
the combined measurements of current and shot noise 
provide a 'spectroscopic' tool to gain information about 
the level structure. For 'local' systems such as single dots 
or two parallel dots (which can be considered as a single 
dot with two or more single-particle levels) the shot noise 
is known to be mostly sub-Poissonian (5 < lei. Super- 
Poissonian noise is predicted in asymmetric situations, 
with level- or spin-dependent couplings to the quantum 
dot la, J4 , | a , or sometimes in the Coulomb blockade re- 
gion 181 |l3l. Similar considerations hold for two dots in 
series 

In this letter we show that starting with a chain of 
three dots the noise may become super-Poissonian even 
if the dot-electrode couplings are fully symmetric (spa- 
tial and spin symmetry). The shot noise is greatly en- 
hanced in absolute magnitude due to a strong compe- 
tition between "slow" and "fast" transport channels, if 



the non-local Coulomb repulsion dominates the inter- 
dot hopping. As there is no simultaneous enhancement 
of the tunneling current, the shot noise becomes super- 
Poissonian. Such behavior can be achieved involving only 
electronic degrees of freedom, and over a large bias region 
above the sequential tunneling threshold. The strong en- 
hancement of the shot noise in a chain of lateral quantum 
dots should allow its direct experimental detection. 

Model and technique- We consider a series of quan- 
tum dots, each with one spin-degenerate level. Including 
hopping between the dots as well as intra-dot and inter- 
dot (nearest neighbor) Coulomb interactions we arrive at 
the Hamiltonian H = Hi, + £Zr + -ffdots + -St.l + -Ht.r 
with 

-ffdots = ^ ein ia - t ^ ( c L c j<? + h - c ) 

n irT rij a ', (1) 

where i = 1 . . . N and r = L, R. Here, Hi, and Hr model 
the non-interacting electrons with density of states p e 
(assumed as a constant) in the left and right electrode 
described by the Fermi operators a kcrr , ako-r ■ The "ar- 
tificial molecule" term £/dots describes the series of dots 
with Fermi operators c\ a , ci a , on-site energies and near- 
est neighbor hopping t. U and U nn are the strength of the 
intra-dot and nearest neighbor inter-dot Coulomb repul- 
sion. Tunneling between the leads and dots is modeled by 
Ht,l and Ht,R- We restrict ourselves to situations where 
the leads couple only to the one adjacent dot. The cou- 
pling strength is characterized by = 27r|i[| 2 p e , where 
t\ are the tunneling matrix elements. 

We are interested in transport through a chain of 
quantum dots, in particular in the current I and the 
(zero-frequency) current noise S. They are related to 
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the current operator / 
—i(e/H) J2ika ( " 



(I R - Z L )/2, with X, 
/Ta^ ;crr Cj (T — /i.c.^ being the current oper- 
ator for electrons tunneling into lead r, by I = {I) and 



S 



dt{8I(t)Si(0) + Si(0)SI(t)) 



(2) 



where 5I(t) = I(t) - (I). 

For the calculation of the current I and current noise S, 
we use the diagrammatic technique on the Keldysh con- 
tour developed in Ref. [Fa and expanded for the descrip- 
tion of the noise in Ref. |fj In this approach the dot array 
( -Helots) is treated exactly as an interacting many-body 
system, but is perturbatively (incoherently) coupled to 
the electrodes. We diagonalize the dot Hamiltonian i/dots 
and compute the transition rates W xx ' (forming a matrix 
W) between two eigenstates \ and x' m nrs t order per- 
turbation theory in the coupling strengths T[. Given the 
transition rates W, the stationary probabilities p st (the 
diagonal elements of the reduced density matrix of the 
dot array) can be found from a master equation. By 
making use of the conservation of probability, we rewrite 
the master equation in the form Wp st = v, where W is 
identical to W but with one (arbitrarily chosen) row xo 
being replaced with (r, T). The vector v is defined as 
v x = TS XXo . The current and shot noise are given by: 



I= 2h e 



wV 



S = ^ e T (W 77 p st + W 7 PW 7 p st ) 



(3) 



(4) 



The vector e is defined as e x — 1 for all x and the matri- 
ces W 7 (W 77 ) denote the transition rates with one (two) 
current vertex replacing one (two) tunneling vertex due 
to i?r,L or Ht,r. The "propagator" P is obtained from 
WP = Q with Q x , x = (pj - S X , X ){1 - cVxJ. Further 
details of the technique can be found in Ref. |(J 

Note that the noise expression Eq. 0] consists of two 
terms: The first term is associated with W 77 , i.e. noise 
diagrams with two current vertices in a single irreducible 
block. This term typically gives a noise contribution of 
el. The second term is due to reducible noise diagrams. 
It involves the "propagator" P of the dot system and 
therefore accounts for the electronic structure and the 
correlations on the central system. It is the second term 
that is responsible for the interesting correlation physics 
dicussed below. 

Results- In the following we discuss current and shot 
noise for a model of type Eqs. [I] with N = 3 dots and a 
half filled ground state (i.e. one electron per dot). We use 
equal couplings, = Tf = 2.5^eV, so T = rf + T§ = 
5/LteV. Our perturbation expansion is valid for tempera- 
tures much larger than the tunnel couplings. Throughout 



this paper, we choose k^T = hT = 0.025meV which cor- 
responds to T « 0.25iG The dot system is characterized 
by the (uniform) level energy e, the intra-dot 'Hubbard' 
repulsion U and the nearest neighbor charge repulsion 
U nn , which we present in units of mcV. Transport is 
achieved by applying a bias voltage Vb, which is dropped 
symmetrically and entirely at the electrode-dot tunnel 
junctions, meaning that the energies of the dot states are 
independent of the applied voltage. The effects of polar- 
ization due to asymmetric or incomplete voltage drops 
are straightforward to anticipate, but would only add 
unnecessary complexity to the results presented below. 
We point out that the three dot system is the simplest 
system that is not pure interface and therefore is the 
minimal model for a truly non-local "artificial molecule" . 
States which are mainly "localized" at the interfacial dots 
are competing in transport with states that have their 
main weight at the middle dot (even though the Hamil- 
tonian remains fully left «-> right symmetric). This com- 
petition can have dramatic effects in the noise character- 
istics that have no equivalent for smaller systems. 
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FIG. 1: Left panel: current /, shot noise 5 and Fano factor 
F vs. bias voltage for a three-dot chain with fcsT = 0.025, 
e = — 10, t — 2, U — 12 and Unn ~ 0.2 (all energies in 
units of meV). The current increases stepwise while the noise 
is sub-Poissonian for the all bias larger than the sequential 
tunneling threshold. Right panel: current I and noise S for 
same parameters as on the left, except that e = —10.8, U nn = 
5. Above Vt ~ 4mV a region of strongly enhanced noise 
appears, whereas the current behaves in a similar manner as 
in the left panel. All current and noise curves are normalized 
to I = (e/2h)T ~ 60pA and So = (e 2 /2h)F ~ 10- 29 A 2 /Hz, 
respectively. 

We first consider a situation where the non-local inter- 
action U nn is small ( U nn = t/10) and obtain the typi- 
cal behavior for a fully symmetric system, see left panel 
of Fig. ^ The current rises (mostly) monotonically in 
steps, the noise also shows steps, but must not increase 
monotonically (thermally broadened peaks around the 
steps are also possible). The Fano factor F = S/2eI 
will fall between values of 1 (Poissonian noise) and 1/2 
for biases larger than all excitation energies (symmetric 
double barrier noise), though in general it will not fall 
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with a monotonous dependence on bias. At small bias, 
eVb <C fcflT, the noise is dominated by thermal noise, 
leading to a divergence Q of the Fano factor. 

If now the non-local Coulomb repulsion U nn is in- 
creased so that U nn > t we observe (right panel of Fig. 
the usual stepwise increase of current and noise except for 
a bias region beyond the sequential tunneling threshold, 
where the noise is strongly enhanced, accompanied with 
a decrease in the current, signaling negative differential 
conductance (NDC). The noise enhancement (and NDC) 
is the result of the occupation of the quadruplet state Q 
(see inset of Fig. 0and discussion below) that, due to the 
nature of its wave function, contributes a " slow" channel 
of transport competing with other " fast" channels. 

It is important to note that the strong enhancement of 
the noise does not require NDC. As an example, in Fig. [21 
we show current, noise and Fano factor for the same pa- 
rameters as for Fig. ^ right panel, except for a different 
e = —10. The current shows generic behavior, i.e. step- 
wise increase and only a tiny NDC around \% = 8.5mV. 
The noise, however, is tremendously enhanced, with the 
Fano factor F > 1 indicating its super-Poissonian nature. 
While F > 1 in itself is not uncommon in the sequential 
tunneling regime 0, 0, El often F > 1 is achieved by 
a suppression of transport in which the current is more 
suppressed than the noise. Here, the current is not sup- 
pressed, but the noise itself is enhanced in absolute mag- 
nitude over a large bias range, before recovering "normal" 
behavior beyond a bias of Vb — lOmV. 
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FIG. 2: Current I and shot noise S vs. voltage for three dots 
k B T = 0.025, t = 2, e = -10, U = 12 , U nn = 5. The noise 
is strongly enhanced in absolute magnitude above Vb ~ 4mV, 
while the current only slightly increases, leading to a Fano 
factor F > 1. This is due to a competition of "fast" and 
"slow" transport channels. The noise scales like (U n „/t) 2 in 
this regime, while the current saturates as t is lowered. Above 
Vb ~ lOmV "normal" behavior resumes, as the "slow" channel 
of transport is "cut short" (see text). 

For such a noise behavior, firstly, strong electron in- 
teractions (U, Unn 3> T) are needed to have the various 
states compete in transport. Secondly, the outcome of 
this competition is determined by the wave functions of 



the competing states that effectively generate state de- 
pendent tunneling transition rates. Finally, the total spin 
of the states in question can differ by more than the elec- 
tron spin 1/2 , so some energetically and spatially possi- 
ble transition rates vanish due to spin selection rules. In 
the present case, the dominance of the non-local interac- 
tion Unn over the hopping t, leads to a strong modifica- 
tion in the spatial distribution of the relevant many-body 
wave functions as compared to the U nn <C t case. Let us 
consider the 9 states with total spin 1 in the q = 2 charge 
sector that split into three triplets (only T\,T^ shown in 
the insert of Fig. [21 If U nn < t, the lowest triplet (Ti) 
will prefer to have electrons on the middle dot, to maxi- 
mize the kinetic energy. On the other hand, if U nn > tit 
prefers to have one electron each on the leftmost and the 
rightmost dot, thus minimizing both intra-dot and inter- 
dot Coulomb repulsion. This change in the nature of 
the lowest triplet wave function is crucial for the " noisy" 
transport. An equally fundamental role plays the total 
spin 3/2 quadruplet, the second excited state in the q = 3 
charge sector. Due to spin selection rules these quadru- 
plet states can only have tunneling transitions to the 
triplet states of the q = 2 (or q = 4) charge sectors. The 
transitions between the lowest triplet and the quadruplet 
form a "slow" channel of transport. The transition rate 
between them is suppressed by a factor ~ (t/U nn ) 2 , as 
the tunneling can happen only at the interface dots, but 
the wave functions of the triplet and quadruplet mainly 
differ by an electron in the middle dot. The current ef- 
fectively alternates between "fast" tunneling sequences 
between the well connected doublet (Di, G) and singlet 
states (Si), and the "slow" sequences between the low- 
est triplet (Ti) and quadruplet state (Q), leading to the 
super-Poissonian noise in the corresponding bias regime. 

In Fig. [2 on the first plateau for 3.3mV < Vb < 4mV 
with sub-Poissonian noise, transport is mainly achieved 
by tunneling events in which the dot system alternates 
between the ground state G with charge q = 3, and the 
lowest singlet Si with charge q — 2. The states Si and 
the lowest triplet Ti (q — 2) are only split by the small 
"exchange energy" of size J ~ t 2 /U ~ 0.33meV and the 
first excited doublet D\ (q — 2) falls energetically in be- 
tween Si and Ti, see the inset in the right panel of Fig. [21 
Therefore, as soon as the sequential tunneling threshold 
V seci = 2(Es 1 —Ec)/e ~ 3.3mV is overcome, the states Ti 
and Di also participate in transport, though their occu- 
pation turns out to be small (less than 1 percent). Since 
V4cq > 2(Eq — ET x )/e even the quadruplet Q is occu- 
pied, though with the same low probability as the triplet 
Ti 0. However, above V b = 2{E Tx - E G )/e = 4mV 
the transport becomes noisy, with a Fano factor F > 1, 
though the current itself still increases by about 2 per- 
cent (not visible in the figure). When direct transitions 
between the ground state G and lowest triplet Ti become 
possible, the occupation in T%, D\ and Q increases at the 
cost of G and Si , until each state has equal probability. 
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Due to their larger spin multiplicity this means that the 
triplet and the quadruplet states compose about 7/12 of 
the total probability. For even larger bias, states in the 
q = 4 charge sector come into play, which leads to several 
smaller features, until at about Vf, = 8mV the quadru- 
plet becomes even more occupied via the lowest triplet in 
the q — 4 charge sector and the noise and the Fano fac- 
tor reach their maximum. The maximum value of noise 
(and Fano factor) behave as (U nn /t) 2 , whereas the cur- 
rent value at the corresponding plateau saturates as the 
hopping t is decreased. This is due to the homogeneity of 
our Hamiltonian: site disorder, i.e. dot-dependent e, will 
lead to a decreasing current with decreasing t ||. The 
size of t, however, is restricted to the condition t ^> T r , 
which is needed for the division of the total system into a 
coherent dot system with perturbative coupling to elec- 
trodes to make sense. 

At Vb = lOmV transitions between the quadruplet and 
the second triplet (T2) with q = 2 become possible. In 
contrast to the lowest triplet Xj , the states of the second 
triplet (as well as the quadruplet states) do not benefit 
from the hopping t, therefore the structure of their wave 
functions is independent of both the interactions U, U nn 
and the hopping t, only the actual energy varies (Et 2 = 
2e + U nn , Eq = 3e + 2U nn ). The second triplet is spread 
uniformly over the three dots, its transition rate with the 
quadruplet is not suppressed with t/U nn . Therefore, the 
slow channel between the T\ triplet and the quadruplet 
is " cut short" , the transport becomes more homogeneous 
and the Fano factor drops below unity. 

We conclude with a few remarks: 1) As the above ef- 
fect relies on spin quantum numbers we believe that it 
will be robust to standard relaxation processes involving 
phonon and photon emission. A strong magnetic field 
will modify the details, but not the generic behavior of 
the transport. 2) With a change of e (gate voltage), i.e. 
e = —11.5, the enhancement of shot noise can be achieved 
even at the first plateau. Although the details of how 
each state exactly contributes to the transport at a given 
bias are changed, the main mechanism (the slow triplet- 
quadruplet channel) to the enhanced noise is the same. 
3) For U nn <ia similar spatial distribution of the lowest 
triplet wave function can be achieved by raising the level 
position €2 of the middle dot more than t above the e of 
the interfacial dots. However, such a level detuning will 
strongly reduce both current and shot noise, similar to 
the "local" models of Refs. 4) In the Coulomb block- 
ade region and close to the sequential tunneling threshold 
we often find super-Poissonian noise. This enhancement 
is due to the thermal occupation and following sequential 
depletion of excited states that lead to small cascades of 
tunneling events interrupted by long (Coulomb) block- 
ages, resulting in a noisy current. Since the origin of this 
effect is completely different to the one discussed above 
we refer the reader to the work of Belzig and co-workers 
H El who have recently discussed this effect for systems 



restricted to a singly occupied ground state in some de- 
tail. The effect can be even more pronounced in more 
complex systems and ground states with higher charge 
q [14j as the number of available excited states typically 
increases with q. 

In summary, we have discussed the interplay of 
Coulomb interactions and non-locality in fully symmet- 
ric models for a chain of quantum dots. We find super- 
Poissonian shot noise above the sequential tunneling 
threshold where the electronic current is not suppressed. 
The Fano factor becomes large in this regime, because 
the shot noise itself is strongly enhanced. This effect is 
caused by a competition of "slow" and "fast" transport 
channels due to the differing spatial structure and total 
spin of the participating many-body states. We expect 
the effect to be robust, thus offering a way to study shot 
noise experimentally in chains of lateral quantum dots in 
semiconducting heterostructures. 
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